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Abstract. Quantum discord is a measure of the quantumness of correlations. After reviewing 
its different versions and properties, we apply it to the questions of quantum information 
processing. First we show that changes in discord in the processed unentangled states indicate 
the need for entanglement in the distributed implementation of quantum gates. On the other 
hand, it was shown that zero system-environment discord is a necessary and sufficient condition 
for applicability of the standard completely positive description of the system's evolution. 
We demonstrate that this result does not translate into useful quantum process tomography. 
Depending on the details of the preparation procedure only absence of any initial correlations 
may guarantees consistency of the process tomography. 



1. Introduction 

What makes correlations quantum? What makes quantum computer tick? There is no definite 
answer to these questions, even if entanglement is a primary suspect. It is crucial in performing 
tasks that cannot be implemented by local operations and classical communications (LOCC), 
its presence in the system is a natural telltale that something non-local is aground, and one 
entangling gate is a part of universal sets of gates [U |2] . 

Nevertheless, the story is not so simple. The "sausage" or "four and a half tatami teahouse" 
states [3] are nine orthogonal pure product states that are shared between the two parties — 
A(lice) and B(ob). Despite absence of entanglement Alice and Bob cannot distinguish between 
them by LOCC. Similarly the DQCl algorithm demonstrates speed up without entanglement 

Quantum discord ^ |6j is one of the measures of quantum correlations. We review its three 
versions and their properties in the next Section. Since its introduction it was found to be 
related to the efficiency of various Maxwell's demons [TJ [8] , quantum phase transitions [9l [10] , 
and quantum state merging |1H I12j . It was also argued that zero discord is a necessary and 
sufficient condition for the evolution of a system to be of a completely positive type \14\ I15j . 
thus being amenable to a relatively simple characterization, in quantum tomography. 

We present two of its uses. First, ability to create entanglement is necessary for execution 
of two-qubit quantum gates even when they are applied to unentangled states and create no 
entanglement. Changes in discord indicate the failure of a local implementation. Second, having 
a zero system-environment discord does not guarantee that the resulting tomographic description 



is an adequate representation of the relevant gate action. Hence zero discord is not directly 
related to a "simple" evolution of the system. 

2. Three types of discord 
2.1. Correlations and entropies 

The Shannon entropy [16j of a classical discrete probability distribution over a random variable A 
with values a, p{a) = pa, is defined by H{A) = — YlaPa^^SPa- Entropy of the joint probability 
distribution over AB, H{AB), is defined analogously. The Bayes theorem relates joint and 
conditional probabilities, 

p{a, h) = p{a\b)p{b) = p{b\a)p{a), (1) 

where p{a\b) is a conditional probability of ^ = a given that B = b. The conditional entropy of 
A given B, 

H{A\B) = Y,PbH{A\b) = - ^p(a, b) logp{a\b) (2) 

b a,b 

is a weighted average of the entropies of A given a particular outcome of B. 

Correlations between two probability distributions are measured by the symmetric mutual 
information. It has two equivalent expressions |16] . 

I{A : B) = H{A) + H{B) - H{A, B) ; J{A : B) = H{A) - H{A\B) = H{B) - H{B\A). (3) 

Quantum-mechanical (von Neumann) entropy [U [21 [T7] of a system with a density operator 
p is defined as S{p) = —iiplogp. It minimizes the Shannon entropy of probability distributions 
that result from rank-1 positive operator- valued measures (POVMs) that are applied to the 
state p. The minimum is reached on a probability distribution Al^ that results from a projective 
measurement 11 = {Ha, a = 1, . . . d}, Yla = ^a^b = ^ab^a, which is constructed from the 
eigenstates of p, 

S{p) = mh,H{A'^), (4) 

i. e., 

S{p) = H{Af), p = Y,PaK, Pa>0, Y.P- = ^- 

a a 

The first expression for mutual information has an obvious quantum generalization, 

I{pab) ■= S{pa) + S{pb) - S{pab), (6) 

and represents the total amount of quantum and classical correlations |18j . 

To obtain a quantum version of J {A : B), it is necessary to determine a conditional state of 
the subsystem B [8]. Given a complete projective measurement 11 on ^4, a quantum definition of 
J follows its interpretation as the information gained about the system B from the measurement 
on A [5], 

j''\pAB)-=S{pB)-S{pB\n^), (7) 

where the conditional entropy is now given by S{pb\^^) ■= J2aPa'^{PB\Ua)- 
The post-measurement state of B that corresponds to the outcome A = a is 

/9B|n, = (Ila (gi iBPABlIa (g) Pa = tv PA'^a-, (8) 

while the state of B remains unchanged, ps = tv aPab = J2aPaPB\na- 



2.2. Discords Di and D2 

Unlike their classical counterparts, the quantum expressions are generally inequivalent and 
I{pab) > {pab) [5^. 'S]. The difference between these two quantities is 

dY^{pab) ■■= S{pa) + S{pB\n^) - S{pab). (9) 

Its dependence on the measurement procedure is removed by minimizing the result over all 
possible sets of 11, resulting in the quantum discord [5] 

DtiPAB) := mm Df-\pAB). (10) 

Similarly, 

JiipAB) ■■= max J"^ (pab). (11) 

It is possible to define the discord when the difference is minimized over ah possible POVM 
[6]. However, unless stated otherwise we restrict ourselves to rank 1 projective measurements. 

An explicit form of a post-measurement state will be useful in the following text. We denote 
this state as p^ = p^ , where the subscript X stands for A, B, or AB, and use the former 
expression if it does not lead to confusion. After a projective measurement H^, the state of the 
system becomes 

p'ab = Y.P-^^^Pb^ (12) 

a 

where pa and p% = PB\Ua given by Eq. ([8]). 

The discord of the state pab is zero if and only if it is a mixture of products of arbitrary 
states of B and projectors on A [5\, 

PAB = ^Pa'[la^ P%, Pa > 0, ^Pa = l- (13) 

a a 

By using this decomposition and properties of the entropy of block-diagonal matrices [20] we 
can identify (pab) = Hp'ab)^ because S{p'j^) = H{A^) and 

S{p'ab) = H{A'^) + S{pb\TI^). (14) 

The discord is not a symmetric quantity: it is possible to have states with Df{pAB) 7^ Df{pAB)- 
Another possibility is to set 

Jf"" := S{pa) + S{pb) - [HiA"^) + S{pB\Ii^)] = S{pa) + S{pb) - S{pUb), (15) 
arriving to the quantum discord as defined in [7] 

D^ipAB) ■■= mm[F(^n) + 5(pB|n^)] - S{pab), (16) 

where the quantity to be optimized is a sum of post-measurement entropies of A and B. Using 
Eq. dH) we see that Di < D2. It is also easy to see that Di = <^ D2 = 0. Using Eqs. (fT2]l and 
()14p we obtain a different expression for D2: 

Df\pAB) = S{pTB)-SipAB). (17) 

Since the definition of the discord(s) involves optimization, the analytic expressions are known 
only in some particular cases. Moreover, typically it is important to know whether the discord 
is zero or not, while the numerical value itself is less significant. 



2.3. Zero discord and 

It follows from Eq. (fT3]) that if the spectrum of a reduced state pA = YlaPaJ^a. is non-degenerate, 
then its eigenbasis gives a unique family of projectors 11 that results in the zero discord for pab- 
Hence, a recipe for testing states for zero discord and for finding the optimal basis is to trace 
out a subsystem that is left alone (i?), to diagonalize pA and to calculate the discord in the 
resulting eigenbasis. 

If the state pA is degenerate, a full diagonalization should be used. For the state with the 
form of Eq. (fTHj) , each of the reduced states p°^ can be diagonalized as 

P% = Y. "-bPa, P'aP'a = P'a- (18) 

b 

The eigendecomposition of the state pab then easily follows. Writing it as 

/0yli? = ^^i'ar5"^„®P„^ (19) 
a,b 

it is immediate to see that its eigenprojectors are given by Ha® Pa- Hence, if p_B has a degenerate 
spectrum, but pAB does not, the structure of its eigenvectors reveals if it is of a zero or nonzero 
discord. Hence we established [8j 

Property 1. The eigenvectors of a zero discord state Df{pAB) = satisfy 

PAB\ah) =rah\ah), \ab) {ab\ = Ua Pji . (20) 

□ 

This consideration leads to the simplest necessary condition for zero discord |19j : 

Property 2. If Df'{pAB) = 0, then [pA •S' ^b,PAb] = 0. Hence a non-zero commutator 

implies Df'{pAB) > 0. □ 
Naturally, if the state has zero discord, and the eigenbasis is only partially degenerated, we 

can use it to reduce the optimization space. On the other hand, the diagonalizing basis H* is 

not necessarily the optimal basis H or H that enters the definition of Di or D2, respectively. 

Consider, for example, a two-qubit state 

PAB = \{tAB + ba\®lB + caA®a%), (21) 

where are Pauli matrices on the relevant spaces, X = A,B, and the constants b and c 
are restricted only by the requirements that pab is a valid density matrix. For this state 
Pb = 1/2 and pA = diag(l + 6,1 — b)/2. After the measurement in the diagonalizing basis 
H^ = ((1 + o"^)/2, (1 — cr^)/2) the conditional state of B becomes 

PBim = 1/2, (22) 

and the conditional entropy is maximal, S{pb\^^) = log 2. 

On the other hand, in the basis H^' = ((1 + cr^)/2, (1 — ct^)/2) the probabilities of the 
outcomes are equal, p+ = p- = 1/2, but the post-measurement states of B are different from 
the maximally mixed one, 

psinj = ^(l±ca"), (23) 

so the entropy 5(pB|n^) > ^(pBln^). 

This discrepancy motivates us to define a new version of the discord [8j : 



D^ipAB) ■■= S{pa) - S{pab) + s{pB\nt), 



(24) 



where is the set of eigenprojectors of pA- For the degenerate case it can be introduced 
using the continuity of entropy in finite-dimensional systems [20j. By applying Eq. to the 
subsystem A we find that simultaneously holds the analog of Eq. ()2.2p . 

J3(pab) = /(Pa|), (25) 

and of Eq. (fTT]) . 

D^{pab)^S{p%)-S{pab)- (26) 
We also arrive to the following ordering of the discord measures: 

Di < < D^. (27) 

There are several important cases when the measures of discord coincide. Most importantly, 
they vanish simultaneously: 

Property 3. Di = ^ D2 = ^ Di = 0. 
The proof follows from Eqs.([li]) and 1^. □. 

For pure states the discord is equal to the degree of entanglement, 

Dfi^AB) = S{(I)a) = E{^ab), i = 1, 2, 3. (28) 

Discord is also independent of the basis of measurement if the state is invariant under local 
rotations [5j. Finally, if ^ is in a maximally mixed state, then = 

These coincidences make it interesting to check when the discords Di and D2 are different. 
By returning to the measurement-dependent versions of the discords, we see that 

Df^'ipAB) = Df^pAB) - MA^) - SipA)]. (29) 

Assume that D2 {pab) reaches the minimum on the set of projectors 11, which are not the 
eigenprojectors of pA- In this case H{A^) — S{pa) > 0, so we can conclude that the strict 
inequality Df < D2 holds, because 

DfipAB) < D^{pab) (30) 

=D^{pab) - [H{A^) - S{pa)] < D^ipAB). 

For example, the state of Eq. ([21]), with b = c = ^, satisfies ^ 0.05, w 0.20 and 
0.21. 

2.4- Some applications 

The discord Df{p) has an operational meaning through quantum state merging [IH I12j . 
Different restrictions on the local Maxwell demons that operate on the subsystems, as compared 
to the power of a global demon, lead to the differences in maximal work extraction that are 
determined by the discords and |8j . Further work on the discord in quantum open 

systems can be found in jTH]. 

Since zero discord is thought to represent absence of quantum correlations, it is interesting 
to investigate the following question. Consider a set of pure orthogonal bipartite states, each 
of which may have a different prior probability, with the ensemble density matrix pab- Does 
the value of D[pab) tell us something about the ability to perfectly distinguish these states by 
local operations and classical communication (LOCC)? The answer is negative [8]: there is no 
relation between D[pab) and local distinguishability, as illustrated in Table 1. We use the nine 
teahouse states 



|0) ® |0± 1)/V2, |2) » |1 ±2)/^/2,|l ±2) |0)/V2, |0 ± l)|2)/\/2, (31) 

of [3], and the result [21j that any two orthogonal (entangled or not) states can be perfectly 
distinguished by LOCC to compile the table. 



Table 1. Local measurability vs. discord 



States 



Discord 



Locally Distiguishable 



9 teahouse states, equal weights 
2 product bi-orthogonal states 
2 entangled orthogonal states 
9 teahouse states, unequal weights 



D^ = D^ = 
D^>0 



yes 
yes 



no 



no 



3. Restricted quantum gates 

It is usually thought that quantum computers are potentially faster then their classical 
counterparts because of the ability to create and use entanglement during the computation 
|22] . Nevertheless, recent results indicate that it is discord rather than entanglement which is 
responsible for the speed-up [H [231 EH [25] . 

Using a simple example we show how changes in discord indicate the need for entanglement 
as a resource even if the processed states are unentangled. More general results may be found 
in [23]. The non-classicality of quantum gates can be examined by the amount of quantum 
resources required for the operation of the gate. Our method is to analyze the resources that 
are needed for a bi-local implementation of a gate. 

A CNOT gate can be implemented bi-locally by Alice and Bob if they share one ebit of 
entanglement per gate use [26]. If they can perform only LOCC they cannot implement this 
gate even on a restricted set of unentagled inputs that are transformed into unentagled outputs. 
One such set L is in Table [2j 



Here cry\Y±) = ztlKj-), ax\X±) = ±\X±), where (Tx,y,z are Pauli matrices. 

We show that the ability to implement the CNOT gate on L without shared entanglement 
allows one to unambiguously discriminate between these non-orthogonal states using just one 
input copy, which is impossible ^17j. Without specifying the local operations of Alice and Bob 
we classify them according to their action on the state 11+). An operation $ is flipping (F) if 
up to a phase $(|y+)) = 1^-), non-flipping (N) if <I>(|y+)) = and undetermined otherwise. 

Knowing the type of the operation allows Alice and Bob to narrow down the list of the 
possible inputs: e.g., Bob's F is incompatible with having input 6, while for Alice's operation 
not to have a definite type excludes both c and d. If one of the operations done is neither F nor 
N, then the type of other operations allows to determine the input uniquely. 

Any pair of outputs can be reset to the original inputs by local unitaries and resent through 
the gate. For example, if the overall operation was of FF type, the operation a^^a^ transforms 
the outputs V'a = 1 1)1^-) and ip'^ = into the inputs V'a and V'c- 

The operations that implement the gate this time may be of the same type as before, or 
different. If the gate is such that there is a finite probability of having a different operation 
type, it will be realized after a finite number of trials. This other type (FN or NF in the 
above example) will uniquely specify the inputs. If the gate's design is such that a particular 
pair of inputs is always processed by the same type of operations, then the gate can be used to 
unambiguously distinguish between one state from this pair and at least one of the two remaining 
states in a single trial. □ 



Table 2. Four inputs/outputs for the CNOT gate 



# State # State 



a |i)|y+)^i|i)|y_) c |y+)|x_)^|y_)|x_) 
b |o)|y+)^|o)|y+) d |y+)|x+)^|y+)|x+) 



The operation in the example changes the discord of the ensemble and has been extended in 
|23j to a more general setting involving a symmetrised type of discord defined as 

D2{p) := m.\n[D^{p),D!^{p)] + 0. (32) 

It is however not yet clear if discord is indeed the best measure for this kind of non locality. It 
is noteworthy to mention that the algorithm described in [4j changes the discord of only one of 
the parties so that D2 remains zero throughout the computation. 

4. (Non) completely positive maps and discord 
J^.l. Completely positive maps and gate tomography 

Environmental interactions are the major obstacle for practical quantum information processing. 
Together with the imperfect tailoring of gate Hamiltonians they are responsible for discrepancies 
between the ideal unitary gates and their experimental realizations. This is why characterization 
of quantum processes is an essential step in implementing quantum technology. The exposition 
below is based on [28]. Experimental realization of this discussion is currently in progress. 

Usually the transformation of input to output states of an open system is described by 
completely positive (CP) maps. Any such map £{p) can be seen as a reduction of a unitary 
evolution of some initially uncorrelated system-environment state pab = PA® where the 
initial states of the system and the environment were pA and w^, respectively. Under the joint 
evolution Uab the state of the system is transformed to 

p'a = £{pa) = tr bUpabU^ = MaPAMl (33) 

a 

where Ma are the Kraus representation matrices [U [2]. Using the spectral decomposition 
CO = '^yPuW){'^\ one recovers the Kraus matrices from 

P = (34) 

This result is based on the absence of prior correlations between the system and the environment. 

The action of a given gate is reconstructed using several methods of quantum process 
tomography. We focus on the standard process tomography, which consists in following. A 
set of input states {pj} is prepared and sent through the process £. From the knowledge of the 
input and the reconstructed outputs £{pk) states the process matrix x P3E71[2U] is reconstructed. 
If Km, m = 1, . . . ,(P form the basis for operators acting on T-Ia, d = dimHAi then 

"^(P) = X] XmnKmpKl. (35) 

m,n 

For a CP evolution all the eigenvalues of x are non-negative, and the Kraus matrices are obtained 
as its (generalized) eigenvectors. 

Various mathematical techniques are used to convert the relative frequencies of different 
measurement outcomes into the states £{pk)- Positivity of states p'j and complete positivity of 
the process £ are enforced, and its violations by the raw data are interpreted as a result of noise. 

The combined state of a system {A) and its environment {B) can be represented in the Fano 
form 

PAB = {^AB + Y ataf CdlB + Y Pj 1a ®'yf + Y1 ^ij'^^ ® ) • (36) 

i j ij 



Here the erf" , i = 1, represent generators of SU(dx), while the real vector a (or a) of size 
— 1 is the generalized Bloch vector of the reduced density operator px- The correlations 
between subsystems A and B are characterized by 

Tij = {jij - ai(3j)/dAdB. (37) 

Presence of correlations not only can lead to a non-CP evolution of the system, but blurs its 
boundary with the environment. Part of the controversy surrounding non-CP maps in literature 
[27j can be traced to this ambiguity as well as to the impossibility of unrestricted creation of 
arbitrary input states of the system (that are to form a tomographically complete basis in the 
standard CP paradigm) pj and their correlations with the environment Tj. 

This situation has an important bearing on the process tomography and action of quantum 
processing devices. In the following we show that while zero discord Da(pab) ensures a 
complete positivity of the resulting evolution [T3], it does not necessarily translate into a useful 
tomographic description of the ensuing evolution. 

4-2. Role of preparations 

State preparation is a vital aspect of quantum process tomography. When initial correlations 
are present in the input state, the preparation procedure used has considerable impact on the 
outcome. We discuss two such schemes: state preparation by a measurement and rotations and 
state preparation by projective measurements |3U| . 

In the former case (Method 1) a single post-selected state that corresponds to a projection 
operator pjj = is used. Once we have performed the projective measurement the required 
input states are obtained by applying the appropriate matrix from a set of unitary rotation 
matrices {Ra} such that RaPHRa = Pa where pa is run over the tomographically complete set 
of input states. 

Hence the preparation maps for Method 1 are given by 

'Paif^Xs) = — {Rallfj 1b) (P^Ab) ^4 ® ^b) =Pa® i^H, (38) 

where pa = tr (H^p'J') is the probability of detecting the outcome corresponding H, and Rj 
brings \H){H\ to one of the states pa = \ipa){ipa\- The post-preparation state of the environment 
is dependent on both the initial state p"^ and the measurement operator Hh- 

State preparation using only measurements (Method 2) utilizes the set of projectors {Ha} 
where H^ = Pa- Hence the preparation procedure for projective measurements is given by the 
collection of maps {Va} where 

Wab) = - ® Ib) f^XB i^a ® 1b) = Pa ^ C^a- (39) 

Pa 

where pa = tr (H^p™) is the probability of detecting the outcome corresponding to Hj. In this 
preparation scheme the state of the environment, ujj depends on the probe state of the system, 
Pk, hence the process matrix will not necessarily be CP. 

Lemma 1: Zero discord of the initial state does not guarantee a CP evolution in a 
tomographic procedure. 

Consider the state 

p% = lm{O\0\O){O\ + \l){l\^\l){l\), (40) 

which satisfies even a stronger requirement D'^{p™^) = D^(p^^^) = 0. If the tomographic 
procedure is based on preparations by measurements, the resulting process matrix x has the 
eigenvalues (1 ± \/3/2, ib\/3/2) and hence the evolution is not CP. □ 



Preparation by measurement and rotations is guarantied to lead to a CP evolution. However, 
we should be concerned that the actual state on which the gate will operate is different from 
the ones prepared in this procedure, and hence the predictive power of the resulting process 
matrix x will be low. Having a zero discord does not alleviate this difficulty, as shown below. 
We consider different reconstruction of process matrices for a generic fixed overall unitary Uab ■ 

Lemma 2: Reconstructed process matrices xVa that follow from the spanning set of 
measurement preparations {V^} that are given by projectors H;^, a = 1, ... — 1, applied 
to the same initial system-environment state {p^Xb)^ equal if and only if this initial state is 
simply separable, (p™^) = (p^) (w^)- 

The sufficient part is obvious. To prove the necessary condition we recall that Kraus matrices 
form an eigendecomposition of the process matrix, hence from Eq. (|34p it follows that all post- 
selected environmental states have to have the same spectrum. Hence up to a local reshuffling 
after the preparation by measurement H^ 

I^XB^r^{f^XB)='^a0U;l (41) 

Using Fano decomposition of Eq. (j36p we find that for all a the post-measurement state of B is 



u% = tr aV^{p%) = ^ (1^ + ^ /3jVf ) + ^ 5^ ( E n^tr (H^a,^) j af , (42) 



hence for any j 



(43) 



for all a in all sets of tomographically complete projectors. As a result, T 
state pXb is indeed a direct product p^ <^ . 



J" = so the initial 

□ 



4-3. Examples 

We illustrate the differences between the outcomes of Methods 1 and 2 by considering qubits 
as the system and environment and a CNOT gate with system A as the target being the 
overall unitary. Let the initial state of joint system AB to be the maximally entangled state 
= where |$+) = (|00) + \ll))/V2. 

The reconstructed process matrices are 



Xi 



( 1 












/ 


2 







1 


\ 














1 










1 


l + i 
















' ^2 ~ 2 




-l + i 


1 


2 







V 1 












v 


1 


1 - i 











(44) 



in the former case we have a 



for Methods 1 and 2, respectively. As follows from Eq. 
CP evolution (actually it is an identity £i = I), while X2 on the other hand has eigenvalues 

(1 + -^,—-^,-^,1 — -^) and hence the reconstructed process is not completely positive. To 
see why this happens consider the state of the environment after each preparation procedure. 
We have that all states are different, 



1 



(45) 



While we used a maximally entangled input state, this result holds even for separable inputs. 
For example, for the input state p^s = ^ (ph ® Pa + Pd ® Pv), the process matrix X2 has the 
eigenvalues {1.642,0.507, —0.253,0.105}, so the evolution is stih non-CP. 



We would like to thank C. Rodriguez-Rosario and K. Modi for their helpful comments and 
suggestions. 
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